Reentrant magnetic ordering and percolation in a spin-crossover system 
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Spin-crossover compounds, which are characterized by magnetic ions showing low-spin and high- 
spin states at thermally accessible energies, are ubiquitous in nature. We here focus on the effect of 
an exchange interaction on the collective properties for the case of non-magnetic low-spin ions, which 
applies to Fe 2+ compounds. Monte Carlo simulations are used to study a three-dimensional spin- 
crossover model for the full parameter range from essentially pure high spin to essentially pure low 
spin. We find that as the low-spin state becomes more favorable, the Curie temperature drops, the 
universality class deviates from the three-dimensional Heisenberg class, and the transition eventually 
changes to first order. A heat-bath algorithm that grows or shrinks low-spin and high-spin domains 
is developed to handle the first-order transition. When the ground state has low spin, a reentrant 
magnetic transition is found in a broad parameter range. We also observe a percolation transition 
of the high spins, which branches off the first-order magnetic transition. 

PACS numbers: 75.10.Hk, 75.20.Ck, 75.50.Xx 



I. INTRODUCTION 

Spin-crossover materials are characterized by ions or 
atoms that can be in either a low-spin (LS) or a high- 
spin (HS) state. If the LS state is the ground state, 
they often show a crossover or phase transition to HS 
at higher temperatures, since the HS state has larger de- 
generacy and is thus entropically favored. 1 Spin-crossover 
systems include organometallic complexes organic 
radicals^ Prussian Blue analogues jLSuLIfi and other in- 
organic transition-metal salts . 11 ' 12 For example, a prin- 
cipal mineral in the Earth's lower mantle is the spin- 
crossover compound Mg^^Fe^O^i 2 - 

If the spins are carried by transition- metal ions, the en- 
ergy difference between HS and LS is due to the interplay 
of crystal-field splitting of the d-levels and Hund's rule 
couplingi 2 -^ For Fe 2+ ions the LS state has spin quan- 
tum number S = 0, whereas the HS state has S = 2. 
Consequently, LS ions are essentially non-magnetic and 
effects due to the LS/HS transition are expected to 
be pronounced. These systems are related to diluted 
spin models ) 13 : 14 ' 15 ' 16 ' 17 but in spin-crossover systems the 
LS /HS degree of freedom is not quenched but dynamical. 

Spin-crossover systems show interesting collective be- 
havior, since they have two coupled degrees of freedom: 
the Heisenberg-type spin orientation and the Ising-type 
LS/HS degree of freedom. In one-dimensional systems, 
quantum fluctuations are important, in particular in the 
case of an antiferromagnetic exchange interaction ! 18 i 19 
At zero temperature quantum effects in this case stabi- 
lize complex phases with magnetic unit-cell lengths of 
3, 5, 7, . . . lattice constants*^ 

In this paper, we study a three-dimensional (3D) spin- 
crossover model including Heisenberg spins with isotropic 
exchange interaction. The LS state has spin 5 = 0. This 
is the 3D version of the spin-crossover chain studied by 
Timm and Schollwock.— We expect quantum fluctua- 
tions to be much weaker in 3D and employ a classical 
description. Nishino and coworkers 2 ^ have studied a 3D 



model with Ising-type exchange, low spin S = 0, and 
high spin S — 1/2, i.e., two possible orientations of the 
high spins, employing mean- field theory and Monte Carlo 
(MC) simulations. Konishi et al— have studied an Ising- 
type model for Prussian Blue analogues. Our model is 
different in that we describe Heisenberg spins, which are 
3D unit vectors interacting via an isotropic exchange 
interaction. 

We perform MC simulations for the Heisenberg spin- 
crossover model on a simple cubic lattice. We are inter- 
ested in the overall phase diagram, not in high-precision 
values for the transition temperature or critical expo- 
nents. In Sec.|TT]we briefly introduce the model and dis- 
cuss the simulation technique. In Sec. lIIII we present and 
discuss our results, which are summarized in Sec. IIVI 



II. MODEL AND SIMULATIONS 

We are interested in 3D systems containing Fe 2+ ions. 
Since the HS state has S = 2 and the system is 3D, we 
expect quantum fluctuations to be relatively weak. We 
therefore start from the classical Hamiltonian 



H : 



-V^CTiO-j 



(ij) 



(1) 



where o~i equals +1 (—1) and |Sj| is zero (unity) for the 
LS (HS) state. We thus absorb the magnitude of the high 
spins into the exchange coupling J. i runs over all sites 
of a simple cubic lattice, while {ij) denotes all nearest- 
neighbor bonds, counting each bond once. We can map 
the antiferromagnetic case J < onto the ferromagnetic 
case J > by flipping all spins on one sublattice. We 
thus restrict ourselves to J > without loss of generality 
and express all energies in units of J. The coupling V 
describes the affinity for equal LS/HS states on neigh- 
boring sites. This interaction is thought to be of mostly 
elastic origin i 22 ' 23 Finally, the Ising "magnetic field" A 
describes the on-site energy difference between HS and 
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LS states. To be precise, the energy difference is 2A and 
A > favors a LS ground state. 

The ground state is a ferromagnetic HS state for 
A/J < 3/2 and V > A/6 - J/2, a pure LS state for 
A/ J > 3/2 and V > -A/6, and a checkerboard LS/HS 
state for V < A/6 - J/2 and V < -A/6. This is sim- 
ilar to the T — phase diagram for one-dimensional 
chains with ferromagnetic Ising exchange interaction^ 
We mostly concentrate on the case V — 0. 

We also need to specify the degeneracies of the LS and 
HS states. In the MC simulation these govern the at- 
tempt rates for LS-to-HS and HS-to-LS flips. Only the 
ratio of degeneracies is important because the physics is 
not affected by an overall constant factor in the partition 
function. For HS S — 2 the underlying quantum system 
has a degeneracy factor of <?hs = 5 per site in the HS state 
and gLS = 1 in the LS state. However, it is well known 
that the ratio G = gas I 3ls can be strongly enhanced 
due to softer vibrations in the HS state ^^^^L 

Our results apply to arbitrary degeneracies since one 
can absorb G into an effective Ising magnetic field— i 28 i 29 

A off = A-^lnG (2) 

(we take ks = !)■ The simulations are therefore per- 
formed for G = 1, i.e., with the same attempt rates for 
LS-to-HS and HS-to-LS flips. For specific systems, A is 
then replaced by A e g. 

We now turn to the MC simulations for this model. We 
simulate finite systems of size L x L x L with periodic 
boundary conditions. The Heisenberg spins Si are stored 
even for the LS sites, but the acceptance probabilities 
do not depend on them. We employ several different 
updates: (a) We use Wolff single-cluster updates^ for the 
Heisenberg spins S.; , but restricted to percolating clusters 
of high spins. This means that we choose a random site 
and, if it is in the HS state, grow a cluster to be flipped as 
in the original algorithm,— but only attempting to add 
neighboring HS sites. These cluster moves mostly avoid 
critical slowing down close to the second-order magnetic- 
ordering transition. 

(b) We employ cluster updates of the LS /HS degree of 
freedom, which are also based on the Wolff algorithm 30 
and use a ghost spin to describe the local effective Ising 
magnetic field A. The approach is similar to RoBler's^ 
Since the Heisenberg spins are not updated, we absorb 
them into effective parameters for the Ising model of the 
(Tj. We first make the LS/HS dependence of the exchange 
term explicit by rewriting it as 

(ij) 

With this notation we can use with unit magnitude 
regardless of whether a site is in the LS or HS state. 
Apart from a term that is independent of the er,;, the 
Hamiltonian now reads 

H = - E Vij a l a - ^2K l a t (4) 

to) « 



with 

% ee v + jSt-Sj, (5) 

A, EE Aeff-^S,- J2 S J> ( fi ) 
NN j of i 

where the last sum is over the nearest neighbors of i. The 
Hamiltonian Q describes an Ising model with nonuni- 
form interaction in a nonuniform field. We rewrite the 
external field as the interaction with a ghost spin a g , 
which couples to all sites i with coupling Vi g . If we choose 
<jg = 1 , we require Vi g — Aj . The Hamiltonian then reads 

H = -^VijViaj, (7) 

where the sum includes the ghost site. In the simulation, 
we chose a random site and try to add all bonds emanat- 
ing from it. The probability for adding a bond from i to 
j is Pij — max(0, 1 — e AB / T ), where AE is the change in 
energy according to Eq. 0. If j is a normal lattice site, 
we really add it, but if it is the ghost site, we reject the 
entire cluster update. We then repeat these steps for all 
newly added sites. If we fail to add more sites but have 
not rejected the update, we flip all Oi on the cluster. This 
algorithm satisfies detailed balance i 31 i 32 

We also include (c) local spin rotations and (d) local 
LS /HS Ising spin flips within the Metropolis algorithm as 
a fallback when the cluster updates are inefficient. For 
a LS site, any spin rotation is accepted, since it does 
not cost any energy. We attempt approximately equal 
numbers of MC sweeps using (a)-(d), where one sweep is 
defined as touching every spin once on average. 

The resulting routine is very robust, but suffers from 
one problem: In a relevant parameter range the mag- 
netic transition is strongly first order. Standard meth- 
ods to deal with first-order transitions such as multi- 
canonical algorithm s 33 ' 34 or entropic samplin g 35 ! 36 fail 
for an interesting reason: The low-temperature phase 
has a high HS fraction and ferromagnetic order, whereas 
the high-temperature phase shows a lower HS fraction 
and no magnetic order. In the high-temperature phase 
the Heisenberg spins are decoupled and typically ran- 
domly oriented. Starting from this phase, local LS/HS 
updates are inefficient because the HS phase is stabilized 
by the exchange interaction, which is zero for isolated 
high spins. However, LS/HS cluster updates also have 
very low probability of reaching the HS phase, since the 
random Heisenberg spins arc energetically very unfavor- 
able in the HS state. Thus we need updates that create 
large clusters of high spins and align them. It is difficult 
to do this efficiently while satisfying detailed balance. 

Our algorithm uses a heat-bath technique that at- 
tempts to insert a layer of high spins which are well 
aligned with their neighbors. In detail, at any temper- 
ature we start by simulating two replicas of the system, 
which only differ in the starting configuration, which is 
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pure HS with perfect ferromagnetic order and pure LS 
with random Sj, respectively. If both phases are (meta-) 
stable and not very close together in configuration space, 
these replicas reach quasi-equilibrium in each phase, but 
do not overcome the barrier between them. We then take 
half of each replica and glue them together, thereby cre- 
ating two phase boundaries. For reasons discussed below 
we cut in planes perpendicular to the (111) direction. 

The resulting system is then equilibrated using the up- 
dates (a)-(d) and one additional type (e) consisting of 
randomly selecting two layers L\ and L 2 that are per- 
pendicular to (111), removing L\, and shifting all layers 
between L\ and L 2 (inclusive) by one primitive lattice 
vector x. We then create new spins in layer L 2 using a 
heat-bath algorithm: At each site i € L 2 , en and Si are 
randomly chosen with the proper Boltzmann probability 
distribution. The reason for choosing planes perpendicu- 
lar to (111) is that each site in L 2 only has neighbors that 
are not in L 2 . Thus we can select each (c,-, S^) indepen- 
dently. If we had chosen layers perpendicular to (100) 
we would have had to create a Boltzmann distribution of 
the entire layer. During the equilibration of the glued- 
together system, the updates (e) grow the more favorable 
phase and eventually remove the less favorable one. 

Finally, we perform measurements for the equilibrated 
system. We wait for at least one MC sweep between mea- 
surements. Since the resulting time series is correlated, 
the errors are estimated using the blocking method*^ 
The most important quantities we measure are the aver- 
age square of the magnetization, 



HS , 



HSj 



(8) 



where the sums are only over HS sites, the average fourth 
power of the magnetization, (M 4 ), and the HS fraction 



-y 
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1 - O r 



(9) 



Note that o % = 1 (-1) for LS (HS). 
The Binder cumulant^ 



Co EE - - - 



3 (M 4 



2 2 (M 2 ) 2 



(10) 



is not as useful for the precision determination of the 
Curie temperature Tc as for the pure Heisenberg model, 
since our model has two length scales, the correlation 
lengths £ of the Heisenberg spins Si and £ CT of the LS/HS 
degree of freedom <7j. We expect £ CT to be a continu- 
ous function of temperature through the magnetic tran- 
sition like the HS fraction 7, as long as the transition 
is of second order. The Binder cumulant can be writ- 
ten as a scaling function C2 = CztyL 1 ' 11 ,L/£ a ), where 
t ee (T — Tc)/Tc and v is the critical exponent of £. 
At T — T c , we have C 2 = C 2 (0, L/^{T C )) so that the 
Binder cumulant does not become independent of system 
size. However, it is still true that the Binder cumulant 



approaches a step function for L — > 00. For our definition 
(fTU)) the step is from unity in the ordered phase to zero 
in the disordered phase. We can thus use C 2 = 1/2 as a 
criterion for Tc that is correct in the limit of large L. 

We also measure the fraction P z of configurations with 
a site-percolating HS cluster that wraps around the peri- 
odic boundaries in the z-direction, where the presence of 
this cluster is determined by a variant of the algorithm 
of Machta et al~2- P z (T) curves for different system sizes 
L are expected to intersect at the percolation transition 
in the limit of large L. Other definitions of the wrapping 
probability, e.g., wrapping in all directions, give the same 
result for the transition^ 



III. RESULTS AND DISCUSSION 

In certain limiting cases, our model becomes either a 
pure Heisenberg model or a pure Ising model. We first 
discuss the Ising limit. A is a constant external field 
coupling to the Ising degree of freedom Oi and if we for 
the moment assume the Heisenberg spins to be frozen, 
they contribute nonuniform magnetic-field and coupling 
terms. The total Ising magnetic field is nonzero unless 
A = J = and breaks the symmetry under <jj — > — Oi 
and thus destroys the Ising critical point. The fact that 
the Ising magnetic field depends on the Si and thus fluc- 
tuates does of course not restore the critical point. 

Now consider the Heisenberg limit. The coupling to 
the LS/HS degree of freedom cr, does not break spin- 
rotation symmetry. If we freeze the Oi we have a diluted 
Heisenberg model, as noted above. The Harris criterion^ 
states that quenched disorder should be irrelevant for the 
critical behavior if the critical exponent v of the corre- 
lation length satifies v > 2/d, where d is the dimen- 
sionality of the system. For the pure Heisenberg model, 
v = 0.7112(5)^ larger than 2/3. Since disorder is thus 
irrelevant, the Heisenberg critical behavior should survive 
for a diluted system if it is not preempted by a first-order 
transition. This is indeed found in Ref. fl7l 

However, the di are dynamical. Thus the Harris cri- 
terion does not apply and different critical behavior is 
possible — we return to this point below. Since spin-ro- 
tation symmetry is not broken explicitly by the coupling 
to Oi , a spontaneous magnetic-ordering transition is still 
possible, though. 

In addition, there can be a site-percolation transition 
of the high spins. Clearly, HS percolation is necessary 
but not sufficient for long-range magnetic order. 

To prepare for the discussion of the phase diagram, we 
plot in Fig. Q]the magnetization squared, (M 2 ), and the 
HS fraction 7 as functions of temperature T for elastic 
interaction V — and various values of A/ J. Recall 
that 2A is the bare energy difference between HS and LS 
states. We find four regimes: 

For A/J < 0.30, the magnetization vanishes at a 
second-order transition. The disordered phase above 
Tc is predominantly HS. A — > — 00 corresponds to the 
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temperature 7/7 



FIG. 1: (Color online) Magnetization squared (solid curves) 
and HS fraction (dashed curves) for elastic interaction V — 
and various values of A/J. For the magnetization curves, A 
increases from right to left. For A/J > 3/2 the magnetization 
vanishes at all temperatures. The results have been obtained 
for system size L = 30, 2000 equilibration sweeps and at least 
10000 measurement sweeps. For several values of A/J the 
data points with error bars are shown. 

pure Heisenberg model, where our results are consis- 
tent with T c /J » 1.457219(4) from high-precision MC 
simulations. 4 - Tc decreases with increasing A. 

For 0.30 < A/J < 0.55, the magnetization vanishes 
at a first-order transition. The disordered phase is still 
mostly HS. 

For 0.55 < A/J < 3/2, the transition is still of first 
order. The disordered phase immediately above Tc is 
now mostly LS, but the HS fraction increases again with 
increasing temperature T > Tc- 

For A/J > 3/2 (exact value), there is no magnetic 
order at any temperature. The system is in a pure LS 
state at T = and the HS fraction increases smoothly 
with increasing temperature. 

Figure [2] shows the phase diagram in the (A, T) plane 
for V — 0. Since we are not interested in high-precision 
determination of transition temperatures, we employ the 
following criteria to map out the transitions: Tc is esti- 
mated from the Binder commulant Ci = 1/2 found by 
bisection for L — 30 with 20000 equilibration sweeps and 
2 17 = 131072 measurement sweeps for each data point. 
(For A/J = 0.25 the resulting Tc is indistinguishable on 
the scale of Fig. [2ja) from the one obtained below from 
finite-size scaling.) The temperature of the percolation 
transition is estimated from the intersection of P z (T) for 
L = 30 and L = 40 found by bisection with 20000 equi- 
libration sweeps and 131072 measurement sweeps. 

The phase diagram Fig. [2]Ja) shows that Tc decreases 
smoothly with increasing A. Note that Tc is already 
significantly reduced for A < 0, where the ground state of 
an uncoupled site is still HS. The low-temperature phase 
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(a) LS/HS splitting A/7 




ferromagnetic HS 



°'0.5 ' 0.55 0.6 
(b) LS/HS splitting A/7 

FIG. 2: (Color online) (a) Phase diagram of the 3D Heisen- 
berg spin-crossover model for elastic interaction V = 0. The 
variables are (half) the bare LS/HS splitting in units of the 
exchange coupling, A/J, and temperature. The solid curve 
with circles shows the magnetic transition. It approaches 
the pure Heisenberg limit (dotted line) for A — > — oo. Open 
(filled) circles represent a second-order (first-order) transition. 
The dotted curve with stars shows the HS percolation tran- 
sition. Also shown are typical lines of constant A for ratios 
G = (7hs/5ls = 1,10,100,1000 of effective degeneracies for 
A/J = 1.75. For larger G, reentrant magnetic transitions 
are obvious, (b) Closeup of the intersection of magnetic and 
percolation transitions in (a). 



remains ferromagnetically ordered for small positive A 
due to the exchange interaction, which favors HS. The 
character of the magnetic transition changes from second 
order to first order at A/J — 0.313(5) and Tc approaches 
zero continuously at A/J = 3/2. 

The HS percolation transition intersects with the mag- 
netic transition where the latter is of first order. It 
then coincides with the first-order transition down to the 
point A/J = 3/2, Tc = 0. The intersection point does 
not have any special properties; the percolation transi- 
tion continues for the undercooled magnetically disor- 
dered state (not shown). It is seen that on approach- 
ing the magnetic transition, the percolation transition 
curves away towards the LS phase. This shows that the 
HS phase is favored even though there is no ferromag- 
netic order yet. This can be attributed to short-range 
ferromagnetic correlations. The fluctuations responsible 
for these correlations are relatively strong although the 
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FIG. 3: (Color online) Magnetization squared (solid curves) 
and HS fraction (dashed curves) for elastic interaction V/ J = 
0.25 and various values of A/J. For the magnetization curves, 
A increases from right to left. For A/J > 3/2 the magneti- 
zation vanishes at all temperatures. The results have been 
obtained for system size L = 30, 2000 equilibration sweeps 
and at least 10000 measurement sweeps. The dotted curve 
shows (M 2 ) for the Heisenberg case. 



magnetic transition is of first order, since the parameters 
are close to the endpoint of the first-order line. Figure 
[2]Jb) shows that the percolation transition is reentrant as 
a function of temperature in the vicinity of A/J = 0.553. 

Recall that arbitrary ratios G — <?Hs/ffLS of effec- 
tive degeneracies can be incorporated by replacing A by 
A eff = A - (T/2) In G. Curves of fixed A in the (A cff , T) 
phase diagram are straight lines of slope —2/ InG. Since 
G > 1 (and often G ^> 1, Ref. |25j), the slope is nega- 
tive. Examples are shown in Fig.[2ja). It is clear that for 
not too small G this leads to reentrant magnetic transi- 
tions: With increasing temperature the system goes from 
paramagnetic (usually LS), to ferromagnetic HS, to para- 
magnetic (usually HS). What happens physically is that 
the increasing temperature entropically favors HS until 
the free energy of the ferromagnetic HS state becomes 
lower than that of the LS state. At higher temperature 
thermal fluctuations eventually destroy magnetic order, 
but HS remains entropically favored. 

We briefly consider the case of V ^ 0. Figure [3] shows 
the magnetization squared and the HS fraction for V / J — 
0.25. Positive V favors neighboring sites both in the LS 
or both in the HS state, thereby stabilizing phases with 
HS fractions 7 « or 7 w 1. We expect this to favor first- 
order magnetic transitions accompanied by a large jump 
in 7, which is indeed seen in Fig. [3l For example, for 
A/J = 0.25 the transition is now of first order, whereas 
it is of second order for V = 0, cf. Fig. [TJ In addition, 
at the same A/J, positive V stabilizes the ferromagnetic 
HS phase, increasing Tc- However, the T = transition 
is at A/J = 3/2 for all V/J > -1/4. 

Figure [4] shows the magnetization squared and the HS 
fraction for V/J = —0.25. Negative V favors neighboring 
sites in different LS and HS states. In fact for V/J < 
— 1/4 the T = phase diagram contains a checkerboard 




temperature TIJ 

FIG. 4: (Color online) Magnetization squared (solid curves) 
and HS fraction (dashed curves) for elastic interaction V/ J = 
—0.25 and various values of A/J. For the magnetization 
curves, A increases from right to left. For A/J > 3/2 the 
magnetization vanishes at all temperatures. This is a special 
point in the T — phase diagram where ferromagnetic HS, 
pure LS, and checkerboard phases meet. The results have 
been obtained for system size L = 30, at least 2000 equilibra- 
tion sweeps and at least 10000 measurement sweeps. For two 
values of A/J the data points with error bars are shown. The 
dotted curve shows (M 2 ) for the Heisenberg case. 



LS/HS phase at A/J around 3/2. At the point V/J = 
— 1/4, A/J = 3/2 three phases (ferromagnetic HS, pure 
LS, and checkerboard) meet. 

We find two phase transitions for intermediate values 
of A/ J: The HS fraction 7 shows a sharp kink above Tc- 
We also note that the HS fraction 7 is close to 1/2 for 
T >Tc, consistent with checkerboard order of the LS/HS 
degree of freedom. An analysis of the average staggered 
Ising order parameter, (\L~ 3 J2i(~^T a i\) ( n °t shown) 
finds long-range checkerboard order in the intervening 
phase. This phase thus appears at finite temperatures 
for V/J = —0.25, although at T = it is absent except 
at A/J = 3/2, as noted. 

In addition, first-order transitions are again favored 
compared to V — 0: For A/J = 0.25 the transition 
is of first order for V/J = —0.25 and 0.25 but not for 
V = 0. Furthermore, Tc for the same A/ J is significantly 
reduced compared to V > 0, since neighboring LS/HS 
pairs are of course unfavorable for magnetic order. 

We next discuss the nature of the magnetic transi- 
tion in some more detail. As noted above, the Harris 
criterion^ is not applicable to our model. It is thus rea- 
sonable to ask whether the dynamical dilution changes 
the universality class where the transition is of second 
order. To answer this question, we perform finite-size 
scaling for the case of A/J = 0.25, for which the tran- 
sition is of second order, but close to the end point of 
the first-order line, see Figs. [1] and [^a). The magneti- 
zation squared close to the transition should scale like 
(M 2 ) L 2I3 / V ~ ^{tL 1 ' 1 ') with t = (T — T c )/T c , where 
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FIG. 5: (Color online) Scaled magnetization squared, 
(M 2 ) L 2f3/ ", vs. scaled temperature, tL 1/v , for V = 0, A/J = 
0.25 for system sizes L = 25, 30, 40, 50. The results have been 
obtained using at least 2000 equilibration sweeps and 262144 
measurement sweeps. Inset: The same quantities with the 
same value of Tc scaled with the exponents (3, v of the pure 
Heisenberg model. 




1.1 1.4 1.5 

temperature 77/ 

FIG. 6: (Color online) Probabilities for finding a spanning 
HS cluster in the ^-direction (upper panels) and HS fraction 
(lower panels) as functions of temperature for (a) A/J — 0.60 
and (b) A/J = 0.75 and system sizes L = 25, 30, 40, 50. The 
horizontal dashed line indicates the critical percolating HS 
fraction in the absence of correlations. The vertical dashed 
lines mark the positions of the percolation transitions. 



$ is a universal function. A least-square fit for sys- 
tem sizes L = 40 and L = 50 yields T c /J = 0.8604(2), 
[3 = 0.26(2), and v — 1.30(8). The errors are the ones 
incurred by the least-square fit, which are much larger 
than the statistical errors of (M 2 ). The scaled magneti- 
zation squared is plotted vs. the scaled temperature for 
sizes L = 25, 30, 40, 50 in Fig. [5] It is apparent that cor- 
rections to scaling are relatively large. Nevertheless, the 
behavior is clearly inconsistent with the 3D Heisenberg 
exponents f3 = 0.3689(3) and v = 0.7112(5)4 2 - The in- 
set in Fig. [5] shows that the scaling fails for Heisenberg 
critical exponents. Since Figs. Q] and [2] do not show any 
qualitative change in the second-order magnetic transi- 
tion as a function of A, we conjecture that the critical 
exponents smoothly approach the 3D Heisenberg limit 
for A — + — oo. 

Finally, we turn to the HS percolation transition. For 
uncorrelated <t,, the transition would take place at a 
HS fraction of 7 = 7 p0 = 0.3116081(21). 43 In our case, 
there are short-range correlations between the <7i, even 
for V = 0, due to the exchange interaction J, which fa- 
vors high spins at nearest-neighbor sites. We have seen 
that these correlations shift the percolation transition to- 
wards larger A. Thus there is no reason to expect the HS 
percolation transition to take place at 7 = 7 p o- Figure [S] 
shows the probability P z of finding a HS cluster spanning 
the system in the z-direction as a function of temperature 
for two values of A and several system sizes. The crossing 
of P z curves indicates the percolation transition^ The 
lower panels show the HS fraction 7 and also the value 
7 p o. The HS fraction at the percolation transition, 7 P , is 
inconsistent with 7 p o: For A/ J = 0.75 we find j p w 0.304 
and for A/J = 0.6, closer to T c , 7 P « 0.295. Thus the 



high spins percolate at a lower HS fraction than they 
would for a random distribution, as expected for positive 
correlations between neighboring high spins. 



IV. SUMMARY 

We have performed MC simulations for a 3D spin- 
crossover model with Heisenberg exchange interaction 
and spin S = in the LS state, as realized in Fe 2+ com- 
pounds. A full range of onsite energy differences 2A be- 
tween HS and LS has been explored, from values strongly 
favoring HS to values strongly favoring LS. We have fo- 
cused on the case of negligible elastic interaction. The 
main results are the following: The Curie temperature 
is significantly reduced already on the HS side (A < 0), 
but initially remains nonzero on the LS side (A > 0) due 
to the exchange interaction J, which favors HS. As A 
increases further, the character of the transition changes 
from second to first order and Tc decreases continuously 
until it reaches zero at a specific value of A. A MC 
method involving growing and shrinking LS and HS do- 
mains has been introduced to tackle the first-order tran- 
sition. 

Where the magnetic transition is of second order, it is 
not in the 3D Heisenberg universality class, at least close 
to the onset of first-order transitions. The system also 
shows a HS site-percolation transition, which intersects 
the first-order magnetic transition. Its position is affected 
by short-range magnetic correlations. They also lead to 
a smaller HS fraction 7 at the percolation transition than 
for a random percolation model. 
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For the realistic case of a LS ground state and an ex- 
cited HS state with much larger degeneracy, we find reen- 
trant magnetic transitions. Increasing temperature en- 
tropically favors HS, which stabilizes ferromagnetic order 
in an intermediate temperature range. 
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